Morse potential derived from first principles 



(N 
O 

Oh: 

< 

(N 



Raimundo N. Costa Filho,^'B Geova Alencar,^'[l| Bo-Sture Skagerstam,'^'[l and Jose S. Andrade Jr.^'H 

^ Departamento de Fisica, Universidade Federal do Ceard, 
Caixa Postal 6030, Campus do Pici, 60455-760 Fortaleza, Ceard, Brazil 
^ Universidade Estadual do Ceard, Faculdade de Educagao, Ciencias e Letras do Sertao Central, 
R. Epitdcio Pessoa, 2554, 63.900-000 Quixadd, Ceard, Brazil. 
^Department of Physics, The Norwegian University of Science and Technology, N-7491 Trondhetm, Norway 
Centre for Advanced Study (CAS), Drammensveien 78, N-0271 Oslo, Norway 

We show that a direct connection can be drawn, based on fundamental quantum principles, 
between the Morse potential, extensively used as an empirical description for the atomic interaction 
in diatomic molecules, and the harmonic potential. This is conceptually achieved here through a 
non-additive translation operator, whose action leads to a perfect equivalence between the quantum 
harmonic oscillator in deformed space and the quantum Morse oscillator in regular space. In this way, 
our theoretical approach provides a distinctive first principle rationale for anharmonicity, therefore 
revealing a possible quantum origin for several related properties as, for example, the dissociation 
energy of diatomic molecules and the deformation of cubic metals. 
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The quantum harmonic oscillator (QHO) is certainly 
one of the most celebrated paradigms in quantum me- 
chanics. Among its several important attributes, the 
QHO can be solved exactly and has been consistently 
used to approximate any potential function when Taylor 
expanded around their minima till second order. More- 
over, the fact that bosons can be conceptually modeled 
in terms of a QHO readily explains its broad applica- 
tion, ranging from fundamental physics, in the descrip- 
tion of a quantized electromagnetic field, to condensed 
matter, for vibrational properties of molecules as well as 
phonons in solids. However, anharmonic potentials are 
very often required to mathematically represent physi- 
cal phenomena. For instance, an adequate description of 
the vibrational modes in diatomic molecules must nec- 
essarily allow for dissociation (i.e., bond breaking) of its 
two bounded atomic nuclei. By increasing the energy of 
the molecule, through heating for example, the system 
starts to vibrate till eventually break down into two non- 
bounded atoms. This is an essential feature of diatomic 
molecules that is not compatible with the QHO model. 

As originally proposed by Phillip M. Morse in 1929 T], 
the so-called Morse potential provides a much better de- 
scription for the potential energy of a diatomic molecule 
than the QHO, being usually written as. 



Vm{x)=D{1 



-ax\2 



(1) 



where x is the distance between atoms, D is the well 
depth related to the molecule dissociation energy, and a 
is an inverse length parameter related to the curvature 
of the potential at the origin. As such, this potential 
has been frequently used as an empirical model for an- 
harmonic interactions in the study of a large variety of 
physical systems and conditions, including the rotating 
vibrational states of diatomic molecules [2| , the adsorp- 
tion of atoms and molecules by solid surfaces Q , and the 



deformation of cubic metals [3]. Figure 1 shows that, ac- 
cording to this potential, the energy difference between 
levels gradually decreases as the level number n increases, 
till it reaches zero, where no more bounded states are al- 
lowed. 

Variants of the Morse potential have also been uti- 
lized to investigate the physical behavior of more complex 
materials. As a typical example, the state-of-the-art for 
atomistic modeling of semiconductor surfaces and inter- 
faces is the Abel-Tersoff potential ^ that has the form 
of a Morse pair potential, but with the bond-strength pa- 
rameter depending on its local environment. Also in the 
investigation of thermal denaturation of double-stranded 
DNA chains, the Morse potential has been successfully 
applied to model hydrogen bonds connecting two bases 
in a pair (rl-fToj. 

It is the purpose of this Letter to show that the Morse 
potential emerges naturally as an effective interaction 
when a particle is subjected to a harmonic potential in 
a contracted space. This physical situation is substanti- 
ated here in terms of the following quantum operator for 



non-additive translation jlll Il2l|: 



(x|L/^(e)|V) = + e(l + .9(72;))), 



(2) 



where there is no restriction on g{-^x). The action of 
Uj{e) on the bra vector (a;| can therefore be expressed 
as {x\U^ie) = (C/^(e)1'|x))t = {x + e{l + g{-fx))\. For 
infinitesimal transformations we obtain that. 



U.y{5x) = 1-1- —p^6x, 



and 



-{x\p^\ij) 



{I + g{"fx)) — {x\ip) , 



(3) 



(4) 



where is the momentum operator. Considering the 
particular case where the function g{'^x) = and a 
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finite displacement a, we can rewrite Eq. ([2]) as, 



{x\U^{a)\^) = ijj xe*" + 



7 



(5) 



from which one can immediately recognize the action of 
the dilation/contraction operator xd/dx. Moreover, as 
defined in Eq. ([2]) , the non-additive operator (a) corre- 
sponds to the infinitesimal generator of the g-exponential 
function ^13i] . 



exp Ja) ^ (1 + (1 - q)a)i/(i-') 



(6) 



with g = 1 — 7. Equation ([5]) represents a fundamental 
mathematical definition for the generalized thermostatis- 
tics of Tsallis and its applications [l4-21|. 

At this point, it is important to state that the momen- 
tum operator is Hermitian with regard to the following 
scalar product: 



dx 



1 + jx 



■)p*{x)(l){x), 



(7) 



where the range of integration shall depend on the spe- 
cific boundary conditions of the system under investiga- 
tion. Equation ([2]) implies that the action of Uj{a) is 
unitary. Indeed, the measure of integration dx/(l + ^x) 
is invariant under the action of the transformation, x 
y — xe"''^ + {e''°- — and the decomposition of the 

unit operator takes the form. 



dx 



1 + 72; 



\x){x\ 



(8) 



The equation of motion for a particle in the x- 
representation of this dilated/contracted space corre- 
sponds to a time-dependent Schrodinger-like equation in 
the form, 



d 



H^{x,t), 



(9) 



where the Hamiltonian operator is H = p^/2m + V{x), 
and the modified momentum operator can be written for 
short as — —ihDj, with D-y = (1 + "fx)d/dx being a 
deformed derivative in space. Equation ^ can then be 
rewritten as. 



d 

ih—'4>(x, t) 



2m 



D^y^{x,t) 



V{x)ij{x,t), (10) 



or, more explicitly, as 



d 

7(1 + '^^)^ a^V^a;, t) + V{x)iJ;{x, t). (11) 

However, it is more convenient to express Eq. (jlip in 
terms of a simple change of variables, as suggested by 



the decomposition of the unit operator Eq. ([8]) . More pre- 
cisely, we can define a variable ry through the differential 
equation drj/dx = 1/(1 + 70;) with boundary condition 
77(0) = 0, whose solution is 



77 



ln(l -t- jx) 
7 



(12) 



From this transformation, it is important to notice 
that the "canonical coordinate" can be written as a; = 
(exp(777) — l)/7- In this way, the finite interval [0, L] for 
X corresponds to [0, L], with L — ln(l -I- '-fL)/j, for the 
variable 77. Equation (jlip rewritten in terms of the new 
variable rj becomes. 



Zfi^0(77, t) = -|-|^</,(77, t) + V,s{7j)<Pi7^, t), (13) 



where 0(77, t) = ^p{x{i]),t) and V^ff(77) = V{x{r])). As- 
suming that. 



(14) 



the transformation leads us back to a more familiar 
version of the time-independent Schrodinger equation, 

E^ii^) = -^^^7j) + V,sivMv), (15) 

If we now consider the problem of a standing wave in 
a null potential, V|,ff(r7) = 0, it follows that the standard 
form of the plane wave solution is recovered in terms 
of the transformed variable 77, namely $(7;) = e^"^^^ . 
In the archetypal case of a harmonic oscillator, V{x) — 
^muj'^x'^, the transformed effective potential becomes. 



2r 



(16) 



where u is the frequency of the oscillator. Strikingly, by 
identifying D = 77ia;^/27^ and a = —7, we conclude that 
Eq. corresponds exactly to the expression ([1]) for 
the Morse potential. To the best of our knowledge, this 
is the first time that a connection based on fundamental 
quantum principles is provided between this potential, 
which has been widely utilized as a consistent descrip- 
tion for the vibrational structure of diatomic molecules, 
and the harmonic potential. A physical interpretation for 
this correspondence can be made in terms of a position- 
dependent (effective) mass induced due to the presence 
of a material body in the system [ll| as, for example, the 
case of electrons propagating through abrupt interfaces 



in semiconductor heterostructures 23-2 



The wave function solution for the quantum Morse os- 
cillator (QMO) has been previously determined as [l[. 



$„(z) = A„z^e-5^L^^(z) 



(17) 
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FIG. 1. The effective potential given by Kff (??) = D{1 - e'^'f 
for 7 = (solid line), 7 = 0.1 (dashed line), and 7 — 0.2 
(dotted line). Here we use uj — 1 and adopt atomic units, 
namely h — m = 1. As shown, the anharmonicity of the 
potential depends on 7. The energy levels for 7 — —0.25 
are shown in the inset. The parameter D corresponds to the 
dissociation energy of the diatomic molecule. 
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2mu}/^^h, s = v/2 



1/2, 



where z 

Ll'{z) = (z-2^eVn!) d" (e-^z"+2'') /dz" is the gener- 
alized Laguerre polynomial [2?|, and the normalization 
constant is given by, 



An — 



T{u-n)' 



(18) 



The energies for the QMO can then be calculated as. 



En 



huj [ n + ■ 



1 



2to 



(19) 



This expression clearly indicates that, in the limit of 
small values of 7, the QHO spectrum is recovered. One 
should also note that the energy levels. 



n+l 



En = fu^ 



^272 



(n + l), (20) 



valid for < n < [rnoj /^'^ — IJ , are no longer equally 
spaced, and approach zero depending on 7 and the quan- 
tum number n. In the main plot of Fig. 1 we show the 
form of the potential ([T6l) for different values of 7, while 
in the inset the energy levels are depicted for 7 = —0.25. 
The potential is symmetric in 7, i.e., ^^(7) = F(— 7). The 
energy difference between levels decreases as the quan- 
tum number increases. 

Finally, we show that an uncertainty relation for the 
QMO can be disclosed in the framework of our theoret- 
ical approach. In order to do this, expected values for x 



and p must be calculated. Considering that the following 
correspondence holds between expected values of a given 
operator O in the x— space and 77— space: 



(21) 



and using a procedure similar to the one adopted in 
Ref. [il], we obtain. 



7ft 



n+- 
muj \ 2 



\ I n,: 



1 - 



2 / I 

n+- 



(22) 
(23) 
(24) 

(25) 



so that the uncertainty relation for the QMO can be writ- 
ten as. 



AxAp ^ h[n+^ 



1 - 



1 

"+2 



(26) 



For the ground state, this expression is identical to the 
one previously obtained in Ref. 11], namely AxAp > 
(fi/2)(l -I- 7(2;)). The uncertainty in the QMO is there- 
fore lower than in the QHO due to the anharmonicity of 
the Morse potential. Note that in both expressions for 
energy (1191) and uncertainty (I26p the parameter 7, which 
defines whether the space is contracting (7 < 0) or di- 
lating (7 > 0), appears squared. Consequently, for the 
harmonic potential in a deformed space, the energy is 
symmetric under contraction or dilation. 

It is interesting to note that the Morse potential can 
be associated to a two-dimensional harmonic oscillator 
in polar coordinates for a; = 1, and considering atomic 
units, ?i = TO = 1. In fact, it is this relationship that 
allows one to solve the QMO using ladder operators [i^ , 
and study this problem through supersymmetry [s^l • As 
a consequence, our approach also provides a conceptual 
basis for the map between a ID-QHO and a 2D-QH0 in 
polar coordinates. This mapping is the result of the new 
commutation relation for x and p generated through the 
translation operator introduced in 11| and used here. 

Previous studies have focused on plausible modifica- 
tions on the position momentum 3l|-|33|, so that a min- 
imum length and momentum could be defined for quan- 
tum theory. In particular, Quesne et al. [s^] have shown 
that, if some special generalized deformed commutation 
relations are employed (e.g., [—e~'^,p] = i[e~'^+l3p'^]), the 
Morse potential can be obtained as an effective potential 
of the theory. In all these studies, however, modified com- 
mutation relations are introduced in an ad hoc. manner, 
i.e., they are not obtained from a first principle mecha- 
nism, like the non-additive translation operator employed 
in this work. Another important point here is that the 
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parameter 7, responsible for the dilation/contraction in 
the translation, corresponds exactly to the minus value 
of the a parameter for the Morse potential. In the par- 
ticular case of the Hydrogen molecule, for example, the 
numerical value of this parameter is 7 = —a = — 1.4 a.u. 

In summary, we have shown that a one-dimensional 
harmonic potential in a space deformed by the action of 
the operator defined as in Eq. ([2]) can be equivalent to 
the Morse potential in a regular space. As the particle 
travels in a different way in the deformed space, it feels 
the harmonic potential as a Morse potential in regular 
space. This equivalence is achieved when 17(72;) = 7X, 
namely for the case in which the translation occurs as a 
contraction in a deformed space, 7 < 0. Such a phys- 
ical framework has perfect analogy with the behavior 
of a position-dependent (effective) mass particle on a 
non-homogeneous substrate, typified by electrons moving 
through abrupt interfaces in semiconductor heterostruc- 
tures. In this particular situation, the anharmonic fea- 
ture of the Morse potential emerges naturally from the 
standard quantum harmonic oscillator. We thus conclude 
that our study, based on a non-additive translation op- 
erator, provides a first principle explanation for anhar- 
monic properties. 
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